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0\ '. Abstract 

-(— > 

s 

'^ ' We present first principles charge- and spin-selfconsistent electronic struc- 

^ , ture computations on the Heusler-type disordered alloys Fes-^Vj^X for three 

Q , different metalloids X=(Si, Ga and Al). In these calculations we use 

the methodology based on the Korringa-Kohn-Rostoker formalism and the 
coherent-potential approximation (KKR-CPA), generalized to treat disorder 
;J] ' in multi-component complex alloys. Exchange correlation effects are incorpo- 



rated within the local spin density (LSD) approximation. Total energy cal- 
culations for Fe3_a;Va;Si show that V substitutes preferentially on the Fe(B) 
site, not on the Fe(A,C) site, in agreement with experiment. Furthermore, 
calculations have been carried out for Fes-^^V^X alloys (with, x = 0.25, 0.50 
and 0.75), together with the end compounds FesX and Fe2VX, and the lim- 
iting cases of a single V impurity in FesX and a single Fe(B) impurity in 
Fe2VX. We delineate clearly how the electronic states and magnetic moments 
at various sites in Fe3_a;Va;X evolve as a function of the V content and the 
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metalloid valence. Notably, the spectrum of Fe3_a:Va;X (X=A1 and Ga) de- 
velops a pseudo-gap for the majority as well as minority spin states around 
the Fermi energy in the V-rich regime which, together with local moments of 
Fe(B) impurities, may play a role in the anomalous behavior of the transport 
properties. The total magnetic moment in Fes^^^V^^Si is found to decrease 
non-linearly, and the Fe(B) moment to increase with increasing x; this is in 
contrast to expectations of the 'local environment' model, which holds that 
the total moment should vary linearly while the Fe(B) moment should remain 
constant. The common-band model which describes the formation of bonding 
and antibonding states with different weights on the different atoms, however, 
provides insight into the electronic structure of this class of compounds. 
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I. INTRODUCTION 

Heusler-type ternary and pseudo-binary compoundgil Y2ZX in the L2i or DO3 structure, 
where Y and Z denote metal atoms and X is a metalloid, display a remarkably rich variety 
of behavior in their electronic, magnetic and transport properties. Among early studies of 
the Fe-based alloys, a note may be made of the work on FeAl,BQ FeSi,lj'cl Fe3Al,Q'Q Fe3Si,all2l 
and some compounds. HjEJ Niculescu et a/.tj give an extensive review of Fes-ajT^jSi alloys 
for various transition metals T. The electronic structure of the Heusler compounds can 
range from metallic to semimetallic or semiconducting. A number of cases of half-metallic 
ferromagnetic phases, where the system is metallic for one spin direction and semiconduct- 
ing for the otherJlj'Ej have been identified. Examples of these are Co2-xFexMnSitIl'E3 and 
Co2MnSii_a;Ge2ji3 and the existence of anti-ferromagnetic ordering in some instances has 
been discussed.lH2l'El It is often possible to substitute on a specific metal site in the lattice 
with other magiietic or non-magnetic atoms, thereby inducing continuous changes in physical 
characteristics. ljEj In view of their tunable magnetic and transport properties, these com- 
pounds have attracted wide attention as potential electronic materials suitable for various 
applications.E3 With all this in mind, it is hardly surprising that Heusler-type compounds 
have been the subject of numerous theoretical and experimental studies over the years. 

Recently, the compound Fe2 VAl has attracted special attention because of the intriguing 
behavior of its electrical resistivity, specific heat, and photoelectric properties.E3 The resis- 
tivity shows semiconductorlike behavior with a negative temperature coefficient suggesting 
an energy gap of ~ 0.1 eV. The photoelectron spectrum on the other hand seems to show 
a Fermi edge which precludes the existence of an energy gap wider than a few hundredths 
of an eV. Finally, low-temperature specific heat measurements for T — > yield a term 7T 
with 7 ~ 14 mJ/mol K^, which results in an effective mass about 20 times as large as 
the bare electron band mass. This mass enhancement is thought to originate from spin 
fluctuationsE3'E3 or from excitonic correlations.E2l It makes Fe2VAl a possible candidate for a 
3d heavy- fermion system. A similar resistivity behavior has been reported for Fe3_2;Va;Ga 
with X near l.Oj^ while Fe2.4V0.6Si shows an onset of the same behavior.EJ 

Concerning relevant theoretical studies, there is a substantial body of literature devoted 
to work on a variety of ordered Heusler-type phases and related systems. Among systems 
of present interest we mention FesSi,!^!^ FesGa,^ Fe3Al,E3E3 and Fe2VAl.EZrEl Little has 
been done on the disordered phases; we are only aware of the study of Fes-^V^jSi alloys in 



Ref. ^, which is based on a non-selfconsistent crystal potential. 

In the present article we report extensive first principles electronic structure compu- 
tations on Fe3_2:Vi,X over the entire composition range for three different metalloids X, 
namely. Si, Ga, and Al. We consider the parent systems FesX for which x = 0, while setting 
X = 1 yields the corresponding compounds Fe2VX. By choosing X=Ga or X=A1 the effect 
of replacing Si by either trivalent Ga or Al has been studied, and by varying x the entire 
composition range of disordered alloys is covered. 

It is well-known that Heusler-type compounds display remarkable 'site-selectivity' 
propertiesEjo in the sense that substituted metal atoms show a preference for entering 
the lattice in specific crystallographic positions. In the generic compound Ya-^^Zj^X, the 
metal atom Z generally prefers B sites in the lattice if Z lies to the left of Y in the periodic 
table, and the A or C positions if Z lies to the right of Y. We have carried out total energy 



calculations which confirm this trend in Fes-a^V^Si. 

An issue of interest concerns the validity of the so-called 'local environmental' model 
(reviewed in Ref. [1^) which holds that the Fe moment scales with the number of Fe atoms 



in the nearest neighbor (nn) shell. The dilution of Fe(B) with non-magnetic V atoms in 
FesSi will then yield a linear decrease in the Fe(A,C) moment, while the Fe(B) moment 
remains unchanged. The reason is that V substitution does not change the nn-shell of 
Fe(B) which continues to contain 8 Fe(A,C) atoms, while the number of Fe atoms in the 
nn-shell of Fe(A,C) decreases progressively. Our computations show substantial deviations 
from this simple picture in Fe3_a;Vj;X, and imply that interactions beyond the nn-shell play a 
significant role in the behavior of electronic structure and magnetic moments. The computed 
spectra also yield insight into a number of other issues, such as the applicability of a rigid- 
band type picture in describing the effects of substitutions in Fe3Si,B and the extent to which 
the ferromagnetism in Fes-ajV^^Xcl can be modeled as a rigid splitting of the paramagnetic 
bands. Finally, we clarify the nature of carriers in Fe3_2;Va;X as a function of composition 
with consequences for transport phenomenaE3'El'E3 in these materials. 

We have attempted to make contact with relevant experiments as far as possible. Al- 
though our primary interest is in the composition dependence of various physical quantities, 
some intercomparisons for the end compounds FcsX and Fe2VX are undertaken. The specific 
experimental data considered are: (i) the composition dependence of the total magnetic mo- 
ment in Fe3_2,.V2,X for the three different metalloids, (ii) the site specific magnetic moments 
in the end compounds FesX and Fe2VX, and (iii) the soft x-ray emission spectra of Jia et 
alx3 on FesSi. Our theoretical predictions concerning the detailed variation of magnetic 
moments on the Fe(B), Fe(A,C), V, and Si sites in Fes-xVx^ show interesting trends which 
should prove worthwhile to investigate experimentally. 

Finally, a few words about our theoretical methods are in order. The disorder is treated 
within the framework of the charge- and spin-selfconsistent KKR-CPA methodology which 
we have developed and implemented in order to handle multi-component random alloys in 
a highly robust manner.lia'EJLJ We use the generalized tetrahedron methodEll to carry out 
k space integrations in disordered muffin-tin alloys, and as a result, our KKR-CPA codes 
allow us to treat the end compounds (FesX and Fe2VX) as well as the properties of single 
impurities in these limiting cases within a consistent, unified theoretical framework. We 
have also generalized the Lloyd formulaPj for the total number of states below any energy 
to multi-component alloys in an analytically satisfactory manner, permitting us an accurate 
evaluation of the Fermi energy in all cases. We are thus able to delineate clearly for the first 
time how the majority and minority spin states and magnetic moments in FesSi develop 
when the Fe(B) position in the lattice is substituted by V atoms, and/or when Si is replaced 
by a metalloid of different valence. The results presented here are highly accurate and involve 
no parameters other than the experimental lattice constants, and constitute a reliable basis 
for testing the underlying KKR-CPA and local spin density (LSD) approximations. 

An outline of this article is as follows. The introductory remarks are followed in Sec- 
tion y by an overview of our KKR-CPA formalism for multi-component alloys. The specific 
formulae used in computing various physical quantities discussed in this article are stated. 
Section |T| outlines the relevant structural aspects of the Fes.^^ V^X compounds. Section [TVI 
summarizes some technical details involved in our computations. Section ^ takes up the 
discussion and presentation of the results, and is divided into a number of subsections in 



view of the quantity and complexity of the material involved. An effort has been made to 
keep the presentation as brief as possible. 

II. AN OVERVIEW OF THE SPIN-DEPENDENT KKR-CPA FORMALISM FOR 
MULTI-COMPONENT COMPLEX ALLOYS 

We consider a multi-component complex alloy where the Bravais lattice is defined by lat- 
tice vectors R„ {n = 1, ..., N} with basis atoms in positions a^ {A; = 1, ..., K}. i^-sublattices 
can be generated from the basis points a^ via lattice translations R„. For simplicity, we 
assume that one of these sublattices, kcp, is occupied randomly by two types of atoms, A 
and B, with concentrations ca and cb, respectively. Other sublattices are taken to be per- 
fectly ordered.E3 The KKR-CPA formalism of interest here proceeds within the framework 
of an effective one electron HamiltonianEa'ErElBEJ where the crystal potential is assumed 
to possess the form of non-overlapping muffin-tin spheres of radius 5*^, i.e., the potential 
is spherically symmetric around each atom and constant (usually defined as the potential 
zero) in the interstitial region. Although such a Hamiltonian is more appropriate for a close 
packed metallic system, a reasonable representation of the crystal potential can often be 
obtained even in open crystals by adding suitably placed 'empty' spheres as basis 'atoms' 
in the lattice. Central for our purpose are matrix elements of the KKR-CPA ensemble- 
averaged Green function G{E), and the one-site restricted Green function G"^^^^ where a 
specific atom Xfc^^ (A or B) sits on the central site of the disordered kcp-th sublattice, 
while all other sites of the knp-th sublattice are occupied by the effective CPA atom. The 
relevant expressions areEJS'&EJ 



< s', r' + a,^,|G^(^)(E)|., r + a^^, >= - E jT\s'r')Z^i^''\sr) + 

aL 
'L' y^^ >^kcpa'L',kcpaL^<^L 



+ E Z,T(-')r.l'f'.,,..l%,r'(-). (2.1a) 



< s',r' + Sikcp\G{E)\s,r + Sikcp >= ca < s' ,r' + Sikcp\G^{E)\s,r + SLk^p > + 



Cb < s ,r +aLkcp\G^{E)\s,r + aLkcp >, (2.1b) 



and 



Wr ' '^7(fc)^ 



<s,r+^,,\G{E)\s,r + e.k>=-E47isr)Z^fi:isr)6,,' + 



aL 



' T ' 



E Z^', ,{sv)TZ,,,^^Z':i{sv), a k and k'^kcp. (2.1c) 



L' ,aL 



Here, r > r, if r > r then J and Z should be transposed. Z^j^ and J^^ are the regular and 
irregular solutions, respectively, of the radial Schrodinger equation within the k-th muffin-tin 
sphere, which may be written compactly as 

E {(E + ^'K's - bf (r)5.'. + t;S'^(r)n • a,,]}zi%r) = 0, (2.2) 



where the up- and down-spin potentials at the fc-th site, v^ and v_ , are combined into a 
scalar part, Vq = (f + +f _ )/2, and a spin-dependent part, t>j = (t>_,_ — 1>_ )/2. n is a unit 
vector along the direction of the magnetic moment, a = {(yxi<^y,(yz) is a vector composed 
of Pauli matrices. Z'^'^^ and J^^"^ are normalized such that outside the muffin-sphere (i.e., for 
r > Sk) they possess the form 

Zfksv) = Y: xAs)ji'{^r)YAr)[r^\^,,^^, - ^^/ExAs)K{VEr)Y,{r), (2.3a) 



45(sr) = xAs)ji{VEr)YL{r). (2.3b) 

Here ^^(f ) is a real spherical harmonic, and L = {I, m) is a composite angular and magnetic 
quantum number index. The spin index a = (+, — ) and the spin variable s = (+, — ) allow 
the treatment of lattices with magnetic ordering. Xcr{s) = 6crs denotes the spin part of 
the wavefunction. h~^{x) = ji{x) + ini{x) is a spherical Hankel function, where ji{x) is a 
spherical Bessel, and n;(x) a spherical Neumann function. The matrix t^''\E) is built from 
on-the-energy-shell elements of the t-matrix of the atom Xk on the k-th site (or A, or B 
atom ii k = kcp)- The elements of t^^\E) are related to the corresponding phase shifts 

n^'hE) by 

^%,.LiE) = -VEexpizrjlS^) sin(r?2V^v'^L'L- (2-4) 

The matrix T"^ (or T^) in Eq. (p.l| ) denotes the so-called central path operator (in the 
sublattice-site representation) for an A or B impurity placed in the KKR-CPA effective 
medium and is related to the medium path operator T'-^^ through the equation 

TMb) ^ j.cP^^ + (^^-1^^ _ r^},)T^r'- (2.5) 

In Eq. ( |2.5D , matrix Tcp is constructed from atomic matrices r^* on sublattices 
with k 7^ kcp, and from the effective scattering matrix t^^ on the fccp-sublattice, i.e.. 

Similarly, for ta[b) we have 

r^ 1 _ f Sk'k [^^'L'L',aL^ {k ^ kcp), , . 

[TAiB^lk'.'L'ML - [ S,, [r^(^)]..,. .„ (k = kcp). ^^-'^ 

The matrix T"-^^ in Eq.( ^.51) is given by the Brillouin zone summation 



T^.'L'ML = ^ E Vcl - B{E, mt'L'MV (2-8) 

with 

[B{E, \.)],,,^,, = Y. exp{tkR^,J[B{E)]^-'^l, (2.9) 

R , 



which are the KKR-complex crystal structure functions]^ defined via a multipole expansion 
of the free electron Green function Go{E): 

< r' + a^. + R„. |G'o(^)|r + a^ + R„ >= 
-i^j:ji{VEr^)ht{VEr^)YL{r')YL{r)6,.,6^,^+j: Y,,{r')[B{E)f-'f^^Y,{r). (2.10) 

L l'L 



The CPA-scattering matrix tcp in the Eq. ( p.5| ) must be obtained by solving the KKR- 
CPA self-consistency condition 

caT^ + cbT^ = T^^. (2.11) 

To solve Eq. ( p.ll|) we use the iteration scheme based on the following expansion 

CB[(Tr)-^ + M-^-(r5p)-T^ (2.12) 



Equation ( p.l2|) allows the computation of T(^p in terms of r^p and T^^ . By carrying out 
the integration of Eq. ( p.8|) , T^ is then determined, and the next iteration cycle can be 



started. This procedure is rigorously convergent and preserves the analytic properties of the 
solutions in the complex energy plancE^TEd This is crucially important because many other 
schemes used in the literature usually fail at some energy points. The problem generally 
becomes more severe as one considers systems with larger number of atoms per unit cell, 
and we have found that Eq. ( p.l2| ) must be the basis of any robust automated procedure for 



obtaining selfconsistent KKR-CPA solutions. 

We have now completely defined the computation of the Green function in Eqs. ( p.l|) 
for a given crystal potential. For carrying out charge- and spin-selfconsistency cycles, one 
other key parameter, namely, the Fermi energy Ep, must be evaluated. In this connection, 
we have developed a powerful version of the Lloyd formulatJa'EJ'El for the total number of 
states below any energy by formally integrating the trace of the KKR-CPA Green function 
exactly in the complex energy plane. The generalization of this formula to spin- dependent 
multi-atom alloys is given below. We emphasize that highly accurate charge- and spin- 
selfconsistent KKR-CPA results of the sort presented in this article would not be possible 
to obtain without the use of this Lloyd-type formula for determining the alloy Fermi energy. 
This is an important point because errors in the Fermi energy determination at any stage of 
the computation impede the convergence of selfconsistency cycles, and degrade the accuracy 
of the final solution for charge- and spin-densities as well as other physical properties. 

We start by taking the trace of the KKR-CPA Green function over the spin- and position 
space, i.e., 

K 

G{E)= Y. Y. d'r<s,v + aik\G{E)\s,v + Siu>. (2.13) 



The integral in Eq. (p. 13 ) extends over the Voronoi polyhedron Vk around the A;-th site. 



and not the muffin-tin sphere, so that space is filled up exactly. Assuming a coUinear 
magnetic structure (same z-axis on each site), the arguments of Ref. |^ can be extended 
straightforwardly to prove that 



G{E)=TrQ + ^{j^ E Tr\n[T^l,-B{EM}- 
-^{ E Trln(0^^) + CATrln(0^) + CBTr/n(0-^)} 

k^kcp 

+ ^{Tr ln[rii - r^^] - c^TrlnirJ^ - r^^] - c^Trln[rji, - r^^]}, (2.14) 

where (j)^i (E) is an energy-dependent renormalization factor for wavefunction Zli \E,r) 
defined by 

zifiE,r) = <l>l^^iE)¥^]\E,r) (2.15) 

with \1/J -^ r' (for r — » 0) . In Eq. ( p.l4| ), TrQ is the free electron contribution, 

TrQ = Y.Y. j d''rji{^r)YL{v)[-t^hi{^r)YL{r)] + ^\n{^fE)K (2.16) 

aL k y^ 

Equation (|2.14]) is cumbersome to use in practical applications as it involves the on-shell 



elements of the t-matrices, r, which do not extend properly into the complex energy plane; 
a form in terms of the logarithmic derivatives turns out to be more useful. First, we write 
the logarithmic derivative at the fc-th mufiin-tin sphere as 



DUE) = Si-\r^Z':,\E,r%=s,- (2.17) 

The diagonal elements of the r-matrices are related to the D-matrices by 



.-m = ^ ' ,, ' ^,^'''^f''^\ (2,18) 

MiVESi.) DUE) - d'^,' MVES,) MVESi,)' 



where 



Di',\E) = Sl^\nUVEr)\^=s,. (2.19) 

We also require the angular-momentum representation of the free-electron Green function 
with position vectors on muffin-tin spheres, 

[Go{E,k)]k'.'L'ML = [D^"^ - D^'^rX'Jk'kh'L + 

j,>iVES,>)[B{E,k)]k'L',kUii^S'')^^'cr (2-20) 

where 

Di?(E) = Sl-^\nhnVEr)\,=s,. (2.21) 

Eliminating r and B in favor of D and G yields our final formula 

GiE) = -^{jj E Tr\n[G^\E,k) + D(^)-Dcp]-'} 

- ^{caTt \n[^^^G^] + cbTt Inf^^^G^] - Tr In G^^} 

+ ^{ E Tr ln[v|>W]}, (2.22) 

k^kcp 



where 

(.AiB) ^ [(^CP)-! ^ ^^^ _ ^^^^^]^ (2.23) 

Equation ( p.22D is formally exact and possesses the form of a perfect derivative. Although 
several terms in Eq. ( |2.22| ) are real on the real axis, their inclusion is crucially important 
for obtaining an analytically correct form which can be used throughout the complex plane. 
Equation ( 2.22|) not only accounts properly for all physical states, but also removes contri- 



butions from spurious singularities present in the scattering matrices and path operators. 
The formal integration of Eq. ( |2.22| ) , 

N(E) = Im / dEG(E), (2.24) 

n J 

— oo 

immediately gives the total number of states, including all core states, below Ep. The 
value of Ep itself corresponds to the condition N{Ef) = Z, where Z is the total number of 
electrons in the Wigner-Seitz cell. 

For a collinear magnetic structure, each spin direction can be treated separately via 
formula ( ^.13| ) yielding the spin-resolved DOS function 



pAE) = ^N^iE). (2.25) 

The magnetic moment /i of a W-S cell can be calculated from 

fi = N+{Ef) - N4Ef). (2.26) 

The spin- dependent charge density at the A;-th site is 

1 ^^ 
pi'\r) = /rfE<a,r + afc|G'(E)|a,r + afc>, (2.27) 

— oo 

where, if A; = kcp, then the k-th atom is taken as an A or B atom and the site-restricted 
Green function of Eq. ( p.laj) is used. The spin density at the k-th site is 

5W(r)=pV')(r)-pL'Hr). (2.28) 

The preceding equations allow the computation of KKR-CPA charge and spin densities 
in the alloy for a starting crystal potential. A new crystal potential may then be constructed 
via the use of the LSD exchange-correlation scheme. The iteration of this procedure leads 
to the fully charge- and spin-selfconsistent KKR-CPA electronic spectrum. 

The magnetic moment on the fc-th site, 

/iW=/iB /'rfVs(^)(r), (2.29) 

is defined as an integral over the muffin-tin sphere volume Qk- Note that within the frame- 
work of the muffin-tin Hamiltonian, this is a unique way of defining site-dependent moments 
in a multi-component system. Since the muffin-tin spheres are not space filling, the sum of 
such individual moments will in general not equal the total moment in the unit cell obtained 
from Eq. ( |2.26|) above, although the differences in the present case turn out to be rather 



small. 
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III. STRUCTURAL ASPECTS 

A brief discussion of the salient features of the crystal structure will help the consideration 
of the electronic properties in the following section. The unit cell shown in Fig. |l| is a cube 
of side a and may be viewed as consisting of four interpenetrating fee lattices denoted by the 
letters A through D; each atom in fact sits at the center of a cube of side a/2 with corners 
occupied by various atoms, so that the packing is identical to that of a simple bcc lattice. 
Fe2VX possesses the classic L2i structure associated usually with the Heusler compounds. 
Here, the two Fe atoms occupy equivalent crystallographic positions A and C, while V sits 
on B sites, and the metalloid X on the D sites. The V atoms on the B sites are surrounded 
by 8 Fe nearest neighbors (nn's) in a bcc arrangement. Each of the Fe atoms in A or C 
positions has 4 V nn's (B) and 4 X nn's (D). The 4 metalloids (D) are located in a relative 
tetrahedral arrangement with respect to each other; this is also the case for the 4 V atoms in 
the B positions. These remarks make it clear that the Heusler compounds contain structural 
units characteristic of metals as well as semiconductors. 

The substitution of V by Fe in Fe2VX giving the disordered alloys Fe3_a;Va;X causes no 
change in the nearest neighbor environment of the V atoms which continue to have 8 Fe 
nn's. The substituted Fe atoms in the B positions, on the other hand, possess 8 Fe nn's as in 
bcc Fe, in sharp contrast to the 4 Fe nn's around each Fe atom in Fe2VX; the end compound 
FeaX is thus rather close to bcc Fe. Although FcsX with two chemical species is classified as 
a DO3 structure, note that in a solid state sense FesX really contains three different types 
of 'atoms' , i.e., two different Fe 'atoms' and the metalloid. It should furthermore be noted 
that, although the alloy Fes-^^VajGa crystallizes in the L2i structure, the end compound 
FeaGa has the DO3 structure only in the narrow temperature range 900 K < T < 920 K. 
Below 900 K the stable phase for small x is the CusAu-type (Llol structure. The DO3 
structure of FesGa is metastable and can be obtained by quenching.Eil'E^ 

The fact that V replaces the Fe atoms only in the B sites of FesX has been adduced from 
NMR, Mossbauer and neutron diffraction measurements (see, e.g., Ref. |T^. More generally. 



in the series Fe3_a;Ta;.Si and Fe3_2.Ta;Ga, where T denotes a transition metal, impurities to 
the left of Fe in the Periodic Table (Mn, V) show a strong preference for the B sites while 
those to the right (Co, Ni) enter at A or C sites. Interestingly, Cr seems to distribute almost 
randomly at A, B, and C sites in Fe3_a;Cra;Si.E3 Reference |3^ has considered the question 
of preferential occupation of various sites in the Fe3Ga matrix via band-theory based total 
energy computations. 

IV. COMPUTATIONAL DETAILS 

We have carried out fully charge and spin self-consistent KKR-CPA computations on 
the series Fe3_2;Va,X, for x = 0.0, 0.1, 0.25, 0.5, 0.75 and 1.0 with the metalloid X 
being Si, Ga, or Al. In the case of the end compounds Fe3X and Fe2VX the KKR-CPA 
results were verified by extensive computations based on our totally independent KKR 
band-structure codes. The self-consistency cycles were repeated for each alloy composition 
until the maximum difference between the input and output muffin-tin potentials was less 
than 1 mRy at any mesh point in the unit cell. Therefore the final potentials used in the 
evaluation of various physical quantities are highly self-consistent. All calculations employ 
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a inaxiniuin angular momentum cut-off lmax='^ and the exchange-correlation functional of 
the Barth-Hedin form.O 

The four basis atoms were placed as follows: Fe(A) = (l/4,l/4,l/4), Fe(C) = (3/4,3/4,3/4), 
Fe(B) or V(Bl=(l/2,l/2,l/2), and X(D) = (0,0,0). The experimental values of the lattice 
const ant scSSe^ of FesX are given in Table |. The lattice constant decreases by less than 
1% in going from FesAl to Fe2VAl.LJ In contrast, the lattice constant slightly increases (less 
than 0.5%) with increasing V content in Fes-^^YxSi. In Fea-^V^Ga also the value of a differs 
by only about 1% between FesGa and Fe2VGa.O Thus, the composition dependence of the 
lattice constant in Fe3_a;V^X is rather weak, and in this work we have neglected the effect 
of this variation, and taken the a value for all compositions to be the same as that of the 
end compound FesX. 

The KKR-CPA cycles were carried out in the complex energy plane using an elliptic 
contour beginning at the bottom of the valence bands, and ending at the Fermi energy 
determined precisely to an accuracy of better than 0.1 mRy via the generalized Lloyd formula 
of Eq. ( p.22| ). This elliptic contour was divided into 12 sections with 4 Gaussian quadrature 
points each, and thus contained a total of 48 energy points; the maximum imaginary part of 



the energy was 0.25 Ry. The KKR-CPA Green function in Eqs. ( p.l|) was computed on a 75 
special k point meshl3 in the irreducible part of the Brillouin zone (BZ) for each of the 48 
aforementioned energy points in order to evaluate the BZ integral of Eq. (|2.8|) . In this way the 
KKR-CPA self-consistency condition was solved at each of the 48 basic energy points to an 
accuracy of about 1 part in 10^, followed by the computation of a new spin-dependent crystal 
potential. The starting potential for the next cycle was typically obtained by a roughly 
10% mixing of the new potential. The solution of the KKR-CPA condition required 1-10 
iterations, while 10-50 charge and spin self-consistency cycles were usually needed depending 
upon the alloy composition to achieve convergence of the crystal potential to an absolute 
accuracy of about 1 niRy. For the final potentials, the total density of states (DOS), site- 
decomposed component densities of states (CDOS), and the /-decomposed partial densities 
of states (PDOS) were computed on a 201 energy point mesh in the alloys, and a 401 
point mesh for the end compounds using a tetrahedral k space integration techniqueEJ (with 
division of 1/48-th of the BZ into 192 small tetrahedra) applicable to the ordered as well as 
the disordered muffin-tin systems. 



V. RESULTS AND DISCUSSION 

A. An Overview within a Simplified Model Density of States 

We present first a relatively simple picture of the component densities of states associated 
with the transition metal atoms in Fes-^VajX. The model density of states of Fig. ^ gives 
the majority (up) and minority (down) spin densities on Fe(A,C), Fe(B), and V sites in the 
limiting cases a; = and x = 1. The positions of the centers of gravity of the Si 3p, Fe 3d and 
V 3d bands are shown, together with the Fermi energies (Ep) for the tetravalent (Si) and 
trivalent (Ga and Al) metalloids. In the following discussion we invoke Fig. |^ frequently in 
order to gain insight into the electronic structure and magnetism of Fe^-xVx^- We emphasize 
however that, even though the model of Fig. ^ captures a good deal of the physical essense 
of the underlying spectrum, the full KKR-CPA selfconsistent results should always be kept 
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in mind. This is especially true in the way Ga and Al-compounds are modeled in Fig. |^ via 
a rigid shift in the position of the Ep compared to the case of Si, the real situation of course 
being more complicated. 

The relative positions of various levels in Fig. ^ can be used to obtain a qualitative 
handle on the movement of bonding and anti-bonding states on various atoms. This aspect 
will play an important role in our analysis below and, therefore, we comment briefly on this 
point with reference to Fig. ^ which describes the so-called common band model of bonding 
of d-hand metals.O Consider two atoms, A and B, with atomic energy level E^ < E'^, which 
are assumed to broaden into rectangular bands of common bandwidth W when the atoms 
are brought together to form a solid. From moment theory, it is known that the center of 
gravity of the local density of states (occupied and unoccupied) must coincide with the local 
on-site energy level Ea{b) (which may be slightly shifted from the corresponding free-atom 
value in order to maintain local charge neutrality). The result is a skewing of the originally 
rectangular local density of states and a new bandwidth Wab- Physically, the skewing 
represents a transfer of charge from the B atom to the A atom until the Fermi levels become 
the same. With this redistribution, the states at the bottom of the band (the bonding 
states) become more concentrated on the A atom while the antibonding states at the top 
of the band are found preferentially on the B atom as shown in Fig. ^. This mechanism 
constitutes a basic ingredient for understanding the electronic structure and magnetism of 
Fe3_i^^V2,.X considered below. A similar discussion based on the bonding of molecular orbitals 
has been given by Williams et alx!\ 

We are now in a position to consider the behavior of the total magnetic moment and 
its constituent parts in Fes-^jV^^X (Figs. H-|]) in terms of the simple model of Fig. 0. Some 
salient features which may be explained are as follows. 

(i) Smaller moment of Fe(A,C) com,pared to Fe(B) in Fe^X. Since the Fe(B) atom is 
surrounded by eight Fe(A,C) atoms in a bcc arrangement, it is not surprising to find the 
component density of states (CDOS) for the d electrons of Fe(B) in Fig. ^ to show the 
familiar structure of two peaks (bonding and anti-bonding) separated by a valley of low 
density of states found in bcc metallic Fe. For the up-spin electrons both peaks are occupied, 
while the exchange splitting pushes the down-spin antibonding peak above the Fermi level, 
resulting in a large magnetic moment on the Fe(B) site. 

Concerning the Fe(A,C) moment, note first that Fe(A,C) is coordinated with four Fe(B) 
and four metalloids, and that the associated CDOS in Fig. ^ possesses extra states between 
the bonding and antibonding peaks. The main difference in relation to Fe(B) is in the 
behavior of the down spins for which Fig. ^ shows that -EFe(A,c) < -SFe(B)- Therefore, 
down-spin bonding states reside preferentially on Fe(A,C) and yield an increased negative 
spin density compared to Fe(B); the corresponding antibonding states on Fe(B) lie above 
Ep and are therefore unoccupied. Fe(A,C) and Fe(B) are quite similar with respect to the 
up spins since bonding as well as anti-bonding states on both Fe sites lie below the Fermi 
energy. The net result is that the total moment on Fe(A,C) is reduced compared to Fe(B). 

(ii) The negative moment of a V-impurity in Fe^X. We see from Fig. ^ that Ey is higher 
than i?Fe(A,c) and its next-nearest neighbor -E'Fe(B) for the up-spin electrons, and thus the 
up-spin bonding states will move away from V sites; the antibonding states will be on V(B), 
but these lie mostly above the Fermi level and are thus unoccupied. On the other hand, for 
down spins Ey lies somewhat below ii^Fe(A,c) and i?Fe(B)- Both effects will tend to induce a 
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negative moment on V impurities. 

(iii) The negative moment of Si in Fe^X. The Si atoms too carry a negative moment, 
albeit small. The reason is that for the down-spin electrons E^i lies well below i?Fe{B), while 
for the up-spin electrons Esi > -EFe(B) • As a result oi p — d hybridization the down-spin 
bonding states will be more predominant on Si while the (empty) antibonding states will be 
found on Fe(B), resulting in a negative Si moment. 

(iv) Positive polarization of Fe(B) impurity in Fe2 VX. Similar arguments explain why 
an Fe(B) impurity in Fe2VSi, at the other end of the concentration range, is positively 
polarized. Figure ^ shows that -E'Fe(B) < -£'Fe(A,c) for the up-spin electrons and thus the 
bonding states will be found predominantly on Fe(B). The opposite situation is true for the 
down-spin electrons. Both effects conspire to produce a strong positive Fe(B) moment. 

(v) Higher moment of Fe(A,C) but not Fe(B) in the Ga- and Al-compound compared to 
the Si-compound. Going to Fe3Ga(Al) we find in Figs. |^ and |^ that, compared to FesSi, the 
Fe(A,C) moments are higher. Figure ^ shows why. Ga and Al are trivalent and therefore 
the Fermi energy is lower. This does not affect the spin-up band which lies in its entirety 
below Ep, nor does it affect the Fe(B) spin-down band because there are hardly any states 
on Fe(B) in the region between the bonding and antibonding peaks. The net effect is a 
reduction of the number of spin-down electrons on Fe(A,C) with a resulting larger moment 
on Fe(A,C), and a negligible effect on the Fe(B) moment. 

The preceding discussion is meant to be illustrative rather than exhaustive. Other fea- 
tures of the behavior of moments in Figs. ^-|^ can be understood at least qualitatively along 
these lines. 



B. Ordered compound FeaSi 

After the introductory discussion of the electronic band structure of the Fe^-xVx^ system 
on the basis of the simple model introduced in the previous section we now turn to the 
calculated density-of-states curves in Fig. |^. A good understanding of FeaSi is essential for 
delineating the effects of substitution on the metalloid and/or the Fe(B) site considered in 
the following sections. Note that the Si 3s and 3p bands show little overlap (Figs. |^c4 and 
1^3) in FesSi, even though in Si the 3s and 3p bands possess a substantial overlap; this is 
the result of an increased Si-Si distance in FesSi compared to Si. 

The bonding between Si and Fe is complex and involves s and p electrons of Si and s,p 
as well as d electrons of Fe. Some manifestations of this bonding are as follows. The Si 
3s states form a semicore band extending below ~ 0.2 Ry; the presence of a finite density 
of states on both Fe sites in this energy region indicates Si-Fe interaction involving Si 3s 
electrons, even though there are no Si atoms in the Fe(B) nn shell. Si-Fe binding via Si 
and Fe p electrons is apparent from the presence of the three-peak structure in the p bands, 
which is most clearly discernable in the Si down-spin p band (Fig. |^3), but is also present 
in the up- and down-spin p bands of both types of Fe sites. Finally, there is the effect of 
hybrid formation between the 3p states of Si and the d states of Fe described by Ho et al\^ 



and discussed in Ref. |33 in connection with FesSi, which tends to concentrate p — d bonding 
states on Fe(B) and enhance the moment on Fe(B). 

The behavior of N{Ef), the density of states at Ep, deserves comment. The states at 
Ep possess mostly d character with some p admixture (Fig. |^). The Fe(A,C) contribution 
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dominates, with the spin-down part being much larger than the spin-up part (Fig. |^al). By- 
contrast, the Fermi level on Fe(B) lies in a fairly low density of states region in the up- as 
well as the down-spin CDOS. The energy dependence of the CDOS in the vicinity of Ep 
is also quite different on various sites. On Fe(A,C), the Fermi level lies near a dip in the 
up-spin CDOS, but in a rapidly decreasing region in the down-spin CDOS (Fig. 0al); a rigid 
upward shift of 0.04 Ry in Ep, for example, would cause a reversal of spin polarization at 
Ep- The situation for Fe(B), on the other hand, is quite the opposite in that a similar shift 
in Ep will induce a rapid increase in the down-spin density (Fig. 0bl). 

Our computed ^-decomposed Si-CDOS (Fig. [7|c) gives insight into the Si L2,3 soft x-ray 
emission (SXE) spectrum of FesSi reported by Jia et a/.E3 The SXE data from FeaSi (see 
Fig. 1 of Jia et al.) display three distinct peaks centered at binding energies of 2 eV, 6 eV, 
and 10 eV. To interpret these results recall first that the L2,3 SXE will only involve s and d 
but not the p partial density due to the A£ = ±1 selection rule for optical transitions. The d 
PDOS (Fig. 0c2) contains many features extending from 0-6 eV below Ep. In Fig. |^ we have 
plotted the sum of s and d PDOS for Si after smoothing the theoretical spectrum to reflect 
experimental broadening.L^ The three peaks in Fig. ^ at binding energies of 1.7 eV, 5.5 eV, 
and 9.2 eV, are seen to be remarkably consistent with the experimental values quoted above. 
In particular, our calculations suggest that the 6 eV peak in the SXE spectrum involves Fe-Si 
bonding d states since both Fe(B) and Fe(A,C) CDOS's possess a substantial density in this 
energy region (Figs. ^1 and_ 
with sp^ bonded Si orbitals. 



djbl, 

|76 



note the scale); in contrast, Ref. ^ associates this peak 



C. Ordered compounds FeaGa and FesAl 

FeaSi and FeaCa are compared first (Fig. ^. The replacement of Si by Ga is seen to 
induce only small changes in the shape of the down-spin CDOS on either Fe(A,C) or Fe(B), 
aside from a relative lowering of the Fermi level to account for the reduced valence of the 
metalloid. The effect on the up-spin CDOS, on the other hand, is more substantial in that 
the shape of both Fe(A,C) and Fe(B) is flatter around Ep in FeaCa compared to FeaSi (e.g.. 
Figs. ^2 andpD2); the Fe-Ga interaction pulls the up-spin states peaking around the Fermi 
level in FeaSi to lower energy. 

The spectra for FeaGa and FeaAl in Fig. |^ are quite similar. As already noted, the Si s 
band in FeaSi (below ~ 0.2 Ry) is more or less core-like; the Ga s band lies at a lower binding 
energy closer to the bottom of the valence band in FeaGa. This progression continues in 
FeaAl where the Al s band overlaps the bottom of the valence band causing a relatively 
greater distortion of states in this energy region (see bottom row in Fig. |^). 

It is noteworthy that substitution on the Si site by Ga and Al not only influences the 
Fe(A,C) but also the Fe(B) CDOS. This is consistent with our observation above in con- 
nection with Fig. 1^ that the metalloid affects the Fe(B) CDOS even though there are no 
metalloid atoms in the Fe(B) nn-shell. It is clear, therefore, that a simple 'environmental' 
type modelEilEHl which is based on taking account of only the composition of the nn-shell 
possesses intrinsic limitations in describing the electronic structure and magnetism of FeaX. 
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D. Ordered compounds Fe2VX (X=Si, Ga, Al) 

The substitution of V for Fe(B) in FeaSi is seen by comparing Figs, ^a and ^a to 
dramatically alter the density of states. The substitution of the high-moment Fe atom by 
nonmagnetic V causes the magnetic moment to nearly collapse. This is reflected in somewhat 
different ways in the up- and down-spin densities. Similar to Fe(B), the V atom in a bcc 
environment displays the familiar d band structure of two peaks separated by a region of low 
density, but since V is nonmagnetic the exchange splitting is close to zero (Fig. p!0|a3). In 
order to accommodate the reduced valence of V compared to Fe, both the up- and down-spin 
antibonding V d states are pushed above Ep. The net changes in relation to FcsSi are most 
clearly visible for the up-spin electrons. The up- and down-spin densities of Fe2VSi possess 
similar shapes with the small magnetic splitting localized essentially on Fe(A,C). The dip 
in the up-spin FeaSi DOS around 0.6 Ry in Fig. |]al, now moved up on the energy scale, 
is partially filled by nnn Fe states around 0.75 Ry (Figs. pUjal and |TD|a2). To the left of 
these there is a complex of Fe-V bonding states followed by metalloid p — d bonding states 
(Figs. p!0|a2 and |l0|a3). Between the Fe states and the antibonding Fe-V states around the 
Fermi level a near gap has formed for both spin directions. 

The character of states at the Fermi level differs greatly between Fe2VSi and FesSi; the 
B site (i.e., V(B) vs Fe(B)) contribution is larger in the former compared to the latter. The 
DOS at Ep is much larger in Fe2VSi and is dominated by Fe(A,C) and V(B) up spins, while 
in FcaSi the Fe(A,C) down spins dominate with other contributions being small. Other 
differences are evident from Figs, ^a and p!0|a; for example, Ep lies in up-spin Fe(A,C) and 
V(B) peaks in Fe2VSi, but in an up-spin Fe(A,C) dip in FesSi. Therefore, we should expect 
Fe2VSi to respond very differently to rigid shifts of the Fermi energy. 

Turning to Fe2VGa (Fig. p!0|b), the spectra are quite similar in shape to Fe2VSi, although 
the minimum in the DOS around 0.8 Ry is somewhat deeper and broader in the Si compound. 
That the spectrum of Fe2VGa does not possess an actual band gap but only a pseudogap for 



either spin direction is seen more clearly from Fig. |TT]. The small downward shift of the Fermi 
level due to the lower valence of Ga places the Fermi level firmly in the (pseudo)gap region, 
thereby precluding moment formation. In sharp contrast to the Si compound, therefore, the 
DOS at Ep in the Ga case is nearly zero for up- as well as down-spins and would increase 
rapidly by a rigid lowering or raising of the Fermi level. 

Notably, the gap between the semicore band around 0.15 Ry and the bottom of the 
valence band in Fe2VSi is larger than in FesSi. This effect is present also in Fe2VGa, and is 
a consequence of changes in the various interactions and not due to a change in the lattice size 
since we have used a lattice constant independent of V concentration in our computations. 
The results for Fe2VAl are seen to be similar to those for Fe2VGa. However, in contrast to 
the case of FesAl (Fig. ^), the bottom of the valence band does not overlap the semicore 
band around 0.25 Ry in Fe2VAl (Fig. [Tgc). 

E. Disordered alloys Fes^V^jX (X=Si, Ga, Al) 
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1. Fe3-^V.,Si 

Fes-ajV^^Si is considered first with the help of Fig. ^. The basic effects outhned in 
the preceding subsection in connection with Fe2VSi are of course at play here. Since the 



driving mechanism is the replacement of Fe(B) by V(B), the b and c panels in Fig. |T2| will 
be considered first. At x = (Fig. |12|bl) the Fe(B) site displays the two-peaked structure 
discussed earlier. A V impurity on that site however shows already the upward shift of the 
two peaks in the up-spin d band. This shift, together with the interaction of the V atom 



with its Fe(A,C) neighbors, has been shown in Subsection |VA| above to lead to the negative 
moment on the V impurity. 

As the V concentration increases, the Fe(B) density of states shows the familiar blurring 
caused by disorder scattering in the alloy. This effect is however highly non-uniform in that 
some states are broadened much more than other states, which reflects large variations in 
the effective disorder parameter as a function of k, E, and spin polarization. lZI Also, with 
increasing x the down-spin bonding peak at 0.7 Ry decreases in size (Fig. |12]b), leaving an 
even larger uncompensated moment in the up-spin band. In the dilute Fe impurity limit 
(Fig. |1^5) the moment on this Fe(B) atom is found to be 3.08 /i^- It is interesting to 
see that in this limit the Fe(B) up-spin bands are virtually undamped again; they have the 
character of relatively sharp impurity levels. 

As a result of the interaction with the V and Fe atoms on the B sites, the Fe(A,C) CDOS 
undergoes quite substantial changes. With increasing V content, the highest occupied peak 
around 0.8 Ry in the up-spin Fe(A,C) CDOS moves to higher energies, which helps deepen 
the low density of states region near the Fermi level and pushes Ep to higher values. These 
large movements in the spectral weights are not present in the down-spin Fe(A,C) CDOS 
although the development of the pseudogap takes place in this case also. The net result of 
the spectral weight shifts is a decrease of the up-spin moment and a simultaneous increase of 
the down-spin moment due to the upward shift of the Fermi level. The two effects conspire 
to cause a rapid depolarization of the Fe(A,C) sites. For x > 0.5 the aforementioned effects 
more or less saturate and there is little further change in the various moments. 

The magnetic moments associated with different sites are seen from Fig. |^ to deviate 
from straight lines joining the x = 0.0 and x = 1.0 values. The moments change essentially 
linearly for x < 0.5 with the Fe(B) component increasing while the Fe(A,C) component and 
the absolute value of the negative moment on V(B) decrease. For larger V concentrations, all 
contributions are nearly flat. Essentially, as Figs. |12| and ^ show, there are three interfering 
mechanisms in going from FeaSi to Fe2VSi: (i) Replacement of high-moment Fe(B) by low- 
(negative)-moment V(B) atoms; (ii) a gradual upward movement of the up-spin Fe(A,C) 
antibonding orbitals which reduces the Fe(A,C) moment, and (iii) in going from FesSi to 
Fe2VSi, as the number of Fe(B) atoms decreases, the Fe(A,C) atoms play an increasingly 
important role in bonding. The concentration dependence of the magnetic moments and 
their lack of linearity results from an interplay of these factors. 

We have made an extensive comparison of our results on Fes.^^VajSi with those of Ku- 
drnovsky et alS3 Despite some overall similarities, our results differ substantially from those 
of Ref. |3^. For example, the Fermi level in our case lies at or near a dip in the majority 



spin-DOS for up to 50% V (Fig. |T2|a), while in Ref. p3 (see their Fig. 8), the majority spin 



N^{Ep) decreases with increasing energy in the 0% and 25% alloy and is essentially at a 
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mininiuin in the 50% V case. The Fermi level in Ref. ^ rises uniformly with increasing V 
content. In sharp contrast, our Ep changes non-linearly, with the Ep values being rather 
close for 0% and 25% V (Figs. p!2|al and 0a2). The site-dependent CDOS's show differences 



as well. The shape of the Fe(B) majority spin-CDOS is seen from Fig. |T2|b to change consid- 
erably in the energy region lying a few eV's below Ep. In Ref. ^ (see their Fig. 10), on the 
other hand, the Fe(B) majority-spin CDOS possesses a roughly composition independent 
three peaked shape. The Fe(B) CDOS at Ep in our computations is quite small for all V 
concentrations, while in Ref. ^ the minority-spin Fe(B) CDOS at Ep is quite large (see 
their Fig. 9). Turning to magnetic moments, we find (Fig.ffl) the Fe(B) moment to increase 



with V content, whereas Ref. |3^ obtains a decreasing Fe(B) moment. Further, our moments 
on all sites vary linearly up to 50% V and remain virtually constant thereafter, and the 
total moment shows a related break in slope around 50% V. The results of Ref. ^ do not 
display these effects clearly as all moments appear to vary roughly linearly up to 75% V. 
The comparisons of this paragraph make it clear that the charge selfconsistency achieved in 
the present work has important consequences for the electronic spectrum; as already noted, 



the results of Ref. 38 are based on a non-selfconsistent crystal potential. 



2. Fes-x Vx Ga and Fe^^x V^Al 



The evolution of the electronic spectrum of Fe3_a;Va;Ga (Figs. |T^) can be understood along 



much the same lines as Fes-^^VajSi, keeping in mind of course the differences in the spectra 



of the end compounds discussed already in Sections |VB| and |VQ above. Interestingly, there 
are differences between the composition dependence of the magnetic moments between the 
two compounds as seen by comparing Figs. | and ||. In Fe3_a;Va;Si, the computed moments 
on all sites remain essentially constant for x >0.5, whereas in Fes-^^VxGa the moments on 
Fe(A,C) and V(B) continue to decrease, and that on Fe(B) continues to increase with x. 
These effects can be related to the behavior of the underlying spectra as follows. Concerning 
Fe(A,C), note that changes in the Fe(A,C) CDOS for either spin direction are quite similar 
in the Ga and the Si compound insofar as the development of the pseudogap around the 
Fermi level and the shifts of spectral weights in the majority spin CDOS (Figs. |T2|a and |1^ 



are concerned. In Fes-ajV^Si, the Fermi level lies in the down-spin pseudogap for x >0.5. 
But, the smaller valence of Ga (compared to Si) causes Ep in Fes-xV^Ga to be relatively 
lower. The minority spin contribution to the Fe(A,C) moment then continues to increase 
for X >0.5 in Ye^_xV ^Ga,, with a concomitant decrease in the total Fe(A,C) moment. 

The composition dependence of the electronic spectrum as well as the magnetic mo- 
ments on various sites in Fe3_a;Va;Al is quite similar to that of Fe3_a;Va;Ga, minor differences 
notwithstanding. While the moments in Fe3_2;Va;Al are presented in Fig. p, the detailed 
results for Fe(A,C), Fe(B), and V(B) CDOS's are not shown in the interest of brevity. 

F. Other aspects 
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1. Rigid hand model vs common-hand m,odel 

Although some features of the Ga and Al alloys can be understood reasonably via a 
rigid band picture (see discussion of Fig. ^ above), our results indicate that a description 
of the electronic spectrum of Fes-^V^^X in terms of any simple rigid-band-type model is 
generally unjustified. The end compounds FesX and Fe2VX possess quite different spectra 
for any given metalloid X, and the shapes of the CDOS's at various sites change with 
Fe/V substitution. Furthermore, the fact that the up- and down-spin DOS's possess very 
different shapes, particularly in the FcsX limit, implies that the ferromagnetic state cannot 
be described properly to be the result of a simple spin splitting in the form of a Stoner shift 
of more or less rigid paramagnetic bands. However, from our discussion it is clear that the 
common-band model constitutes a fertile framework for obtaining insight in the electronic 
structure and magnetism of these materials. While the exchange splitting on the Fe(B) atom 
constitutes the driving force for magnetism in these compounds, the hybridization between 
different states and the resulting differences between the weights of various bonding and 
antibonding states on different atoms lead to a rich variety of behaviors. 

2. Site selectivity 



As pointed out in Section |T|, the issue of site selectivity in the Heusler-type alloys 
has been the subject of numerous studies ,li3l22rE3ll3 with the attention being focused mostly 
on determining whether other metal atoms when substituted for Fe preferentially occupy 
Fe(B) or Fe(A,C) sites in the lattice. The high electronegativity of Si implies that Si attracts 
electrons from the surrounding Fe atoms. Our calculations indicate that in FcsSi the number 
of electrons inside the muffin-tin spheres of Fe(A,C) is 25.115 against 24.975 for Fe(B), i.e., 
Fe(A,C) averaged over both spin directions is more electronegative than Fe(B). The same is 
true in FesGa and FesAl, even though Ga and Al are less electronegative than both Fe(A,C) 
and Fe(B). The Coulomb energy of the crystal lattice will be reduced if electronegative 
Fe(A,C) is replaced by an element more electronegative than Fe (i.e., Co or Ni) or if Fe(B) 
is replaced by a less electronegative element (i.e., Ti, V, Cr, or Mn), thus explaining the 
observed site selectivity. A more quantitative demonstration of these effects, of course, must 
be based on total energy calculations where the LDA has proved a sound basis for metals,EHl 
compounds, e3'EJ and binary alloys.Ej Accordingly, we have computed the total energy of 
the alloy Fe3_j;Va;Si for a number of concentrations x putting V atoms first on the A or C 
sublattice and subsequently on the B sublattice. Figure |1^ shows the difference between the 
total energy for each of these two situations and the sum of the energies of the constituent 
atoms for x < 0.08. It is clear that the total energy is lowered when the V atom occupies a 
B site. 



3. Electrical resistivity 

Figures 1^14 imply that the type, spin, and number of carriers available for transport in 
the Fe3_a;Va;X alloys depend strongly on the V content as well as on the metalloid valence. 



Focusing on Fea-ajV^^Si first, we see from Fig. |T2| that in FeaSi (topmost row) the density 



of states at the Fermi level is dominated by the down-spin Fe(A,C)-CDOS. Therefore, the 
current in FeaSi will be carried primarily by down-spin electrons associated with Fe(A,C) 
sites with relatively little contribution from Fe(B) and V(B) sites. With increasing V content, 
the Fermi level in Fe3_a;Va;Si moves into the pseudogap in the down-spin Fe(A,C)-CDOS; 
for x=0.50, these electrons are seen to be essentially frozen out of the transport processes. 
In fact, for a;=0.50 few carriers of either spin are available and one expects the material to 
possess a high resistivity. For a;=0.75, the up-spin Fe(A,C)-CDOS is larger, and the up- 
spin V-CDOS begins to increase. In the limiting case of Fe2VSi (Fig. |1^, bottom row), we 
see that the current will be carried mainly by up-spin Fe(A,C) and V(B) electrons. Thus, 
in going from FesSi to Fe2VSi, the carriers change from being dominated by down-spin 
Fe(A,C) electrons to up-spin Fe(A,C) and V(B) electrons, and the material goes through 
a high-resistivity range for intermediate compositions. Experiment confirms this picture: 
measurements of the residual resistivity as a function of x in Fea-^V^^Si by Nishino et alB^ 
show a pronounced maximum at x ~ 0.4. 



The situation with Fe3_a;Va;Ga is seen from Fig. |1J to be similar to that of Fea-^^VajSi, 
except that the compound does not display an intermediate range with few carriers around 
x=0.50. Instead, the carriers continue to be dominated by down-spin Fe(A,C) electrons 
for all compositions, with the total number of available carriers going nearly to zero in the 
x=1.0 limit. There seems to be a curious effect in the x = 0.75 case (Fig. |l^b4) in that the 
Fe(B) CDOS displays a substantial up-spin component, suggesting that this alloy may show 
unusual transport phenomena. 

Returning to the compounds Fe2VAl and Fe2VGa, we note that the Fermi surface of 
Fe2 VGa is very similar to that of Fe2VAl.E3'E3 It consists of three small hole pockets at F and 
an electron pocket at X (see Fig. |Tl|), resulting in a very small number of carriers.E3 Recent 
calculations for Fe2VAl which include spin-orbit couplingpl find a further reduction in carrier 
density. This would make these compounds semimetals or even semiconductors, although 
the latter would disagree with a Fermi cut-off reported by Nishino et a/.c3 in Fe2VAl using 
high-resolution photoemission. Another cause for the high and strongly temperature and 
composition dependent resistivity could be strong electron scattering by spin fluctuations. 
Although we find Fe2VAl and Fe2VGa to be virtually non-magnetic in the V-rich regime, 
we have seen that they may nevertheless contain sizable, relatively isolated, local moments 
in the form of Fe atoms on V(B) positions, accompanied by non-stoichiometry or frozen-in 
antisite defects. The band structure of these compounds consists of a set of rather fiat Fe 
and V d bands closely below and above the Fermi level. A similar situation involving / 
bands is known to lead to heavy-fermion behavior in a variety of intermetallic compounds. 
Also, the temperature dependence of the resistivity in Fe2VAl is qualitatively very similar 
to that of well-known heavy-fermion compounds such as CeCu2Si2 and CeCu6.E3 Thus, the 
resistivity behavior in Fe2VAl and Fe2VGa may be a signature of heavy-fermion behavior in 
a 3d intermetallic compound. It is also interestingto note that the concentration x = 0.6 at 
which Fe3_a; Va;Si shows similar resistivity behaviort^ lies in the concentration range where the 
Fermi level is at or very near the strong dip in the density of states which for x ^ 1.0 develops 
into the pseudogap. This suggests that a low carrier density, possibly in combination with 
strong local moments, is a necessary ingredient for semiconductorlike resistivity behavior 
in these alloys and compounds. In this respect resistivity measurements in Fes-^^VxSi for 
X > 0.6 would be of interest. 
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4- Magnetic moments-comparison with experiment 

We comment now on the magnetic moment data given in Table || for the end 
conipounds.E^t3t3C3L3'Ol3i5E3E3 The experimental data in this table is not meant to be 
exhaustive; we focused on sources which report site-dependent magnetic moment informa- 
tion. FesSi has been investigated most extensively. The situation with the Fe(B) moment 
in FesSi appears to be good in that the experimental values lie around 2.2-2.4 /i^, and are 
in reasonable accord with the theoretical value of about 2.5fiB', the value of l.QfiB is based 
on a tight-binding model computation. A similar level of agreement is also seen generally 
on the Fe(B) site in FesGa and FcsAl. In the case of the Fe(A,C) moment, on the other 
hand, the situation is somewhat less satisfactory as the experimental as well as theoreti- 
cal values are more scattered. Theoretically, the metalloid is predicted to exhibit a weak 
negative polarization, but no experimental results are available. Lack of experimental data 
is also evident for the magnetic moments in Fe2VX; notably, experiments of Ref. ^ give 
a near zero Fe(A,C) moment in Fe2VSi in disagreement with computed values of 0.43 fiB 
(this work), and 0.37 /i^ by Ref. |3^. Interestingly, Endo et alM report antiferromagnetic 
ordering in Fe2VSi; on this basis, an Fe(A,C) moment of 0.5-1.0 /is has been deduced.E3 Of 
course, it should be remembered that our calculations do not include an antiferromagnetic 
ground state. 

Also for the alloys Fe3_j.Va;X we are not aware of any experimental work addressing 
contributions of individual sites to the total moment. In Fe3_a;V^Si (Fig. ^), the computed 
total moments are rather close to the measurements of Ref. 0, although further experimental 
work in the V-rich alloys, where the total moment curve is predicted theoretically to display 
a significant change in slope, would be interesting. The experimental data in the Ga£3 and 
A1e3 compounds (Figs. |^ and ^ extend over a wider composition range and here again 
there is a reasonable level of accord with the theory, even though the experimental points 
are generally lower, especially for the 50% V alloy. We should keep in mind however that 
our calculations assume a perfectly random occupation of Fe(B) sites by V atoms while 
the physical system may display short range ordering and clustering effects in addition to 
possible uncertainties associated with site occupancies. 

VI. SUMMARY AND CONCLUSIONS 

We have discussed the electronic structure of Fe3_2;Va;X Heusler alloys, where the met- 
alloid X is Si, Ga, or Al, on the basis of our all-electron, charge- and spin-selfconsistent 
KKR-CPA computations. Specific compositions calculated include the concentrated alloys 
with 25, 50, and 75 atomic percent V substituted randomly in the Fe(B) site, the end com- 
pounds FesX and Fe2VX, and the limiting cases of a single V impurity in FesX and a single 
Fe(B) impurity in Fe2VX. All calculations were carried to a high degree of selfconsistency 
and are parameter-free excepting the use of experimental lattice constants. The Fermi en- 
ergy was evaluated in all cases with a high degree of precision by using a generalized Lloyd 
formula for mult i- component systems. The use of a tetrahedron-type k space integration 
method, which we have extended to handle disordered muffin-tin systems, allows us to treat 
the ordered end compounds as well as the single impurity limits on an equal footing with 
the concentrated alloys. Our results thus provide a reliable basis for testing the underlying 
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framework of the KKR-CPA and LSD approximations. To our knowledge, little theoreti- 
cal work exists in the literature concerning the electronic structure of disordered phases of 
Heusler alloys, although Ref. ^ has previously considered aspects of Fes-^^VajSi based on 
non-selfconsistent computations. 

Highlights of our results are as follows. An examination of the spin-dependent component 
densities of states (CDOS's) at various sites in Fes-j^V^^Si shows that the replacement of 
Fe(B) by V induces substantial shifts in spectral weights, especially in the up-spin Fe(A,C) 
CDOS; the down-spin Fe(A,C) CDOS and the CDOS for either spin direction at other sites 
- Fe(B) or Si - suffer relatively lesser changes. A pseudogap develops around the Fermi 
energy in the up- as well as the down-spin density of states in the V-rich regime. The effects 
of V substitution in Fe3_a;Vx.Ga are similar to those in Fes-ajV^^Si, although a number of 
differences arise from the reduced metalloid valence and the changes in the metal-metalloid 
interaction. Fes-j^VjjAl, on the other hand, possesses a spectrum rather close to that of 
Fe3_3;Va,Ga, some differences in detail notwithstanding. 

The complex spectral changes in the electronic spectra of Fe3_a;Va;X induced by Fe(B)/V 
substitution and by the effects of metal-metalloid interaction cannot be described within a 
rigid band picture based on one of the end compounds. In a similar way, the ferromagnetism 
of these alloys cannot be viewed as a rigid splitting of the related paramagnetic spectra. 
However, the common band model in which the bonding between atoms is described by 
the formation of a common energy band containing bonding and antibonding states with 
different weights on the participating atoms provides insight into the electronic structure of 
these compounds. 

We have analyzed the composition dependence of the magnetic moments associated 
with individual sites in Fe^^x^x^ in detail, and correlated these changes with those in the 
underlying spectra. In Fe3_a;Vj;Si, all moments vary essentially linearly for x < 0.50, and 
remain virtually constant thereafter. V substitution (up to 50% V) decreases the Fe(A,C) 
moment, but increases the Fe(B) moment; Si possesses a very small (-0.08 /i^) negative 
moment. The V impurity in FesSi possesses a negative moment of -0.62 fiB- The total 
moment in Fe3_a;Vj;Si varies non-linearly over < x < 1.0 with a change in slope at a; = 0.50. 
This behavior differs sharply from that expected on the basis of the 'local environment' 
picture which assumes the Fe(A,C) moment to decrease linearly over < s < 1.0 and the 
Fe(B) moment to remain constant. On the other hand, the moments in Fe^-xYxGa and 
Fe^-xVx-^^ show a different type of non-linearity in that the Fe(A,C) moment decreases 
rapidly between 0.5 < x < 0.75, and the aforementioned 'saturation' effect in Fe3_a;Va;Si for 
X > 0.50 is not seen. These results show clearly that the metal-metalloid interaction has 
an important bearing on the magnetism of Heusler alloys. Some measurements of the total 
magnetic moment in Fe3_3;Va:Si, Fe3_a;V^Ga, and Fes-^^V^Al are available, and in this regard 
our theoretical results are in reasonable agreement with the experimental values. Further 
experimental work particularly to obtain site decomposed moments in Fe3_a;Va,X should 
prove worthwhile. Finally, we have delineated the nature of carriers involved in transport 
processes in Fes^^^x^- In the a; = limit, most of the carriers originate in the down-spin 
Fe(A,C) sites. With increasing V concentration, as a pseudogap forms around the Fermi 
energy, the total number of carriers in Fe3_2:VxSi decreases rapidly. In the V-rich alloy, 
the number of carriers begins to increase once again, but it is now dominated by up-spin 
electrons. By contrast, the V-rich compounds in the case of Fe3_a,V2;Ga or Fes-^jV^^Al are 
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predicted to be semimetallic, and to be dominated by down-spin carriers throughout the 
composition range. These dramatic changes in the number, spin, and type of carriers in 
Fes-a^Vj^X and their interaction with strong local moments of Fe(B) impurities may play an 
important role in the anomalous behavior of the resistivity and other transport properties 
and in the possible heavy-fermion character of some of these compounds. 
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FIGURES 

FIG. 1. Crystal structure of Fes-^V^X. The four crystallographic positions, denoted by A-D, 
are shown by shadings of different intensity. The A and C sites are equivalent in the present case. 
Fe and V atoms occupy B sites, while Si, Ga, or Al atoms sit at D sites. 

FIG. 2. Schematic diagram of the up- and-down spin component densities of states on Fe(A,C), 
Fe(B), and V(B) sites for a) FeaSi and b) Fe2VSi. The horizontal dashed lines mark the Fermi 
levels for FesSi and Fe2VSi, the chain lines for the corresponding Ga and Al compounds. The kinks 
in the latter represent an overall shift of the entire set of down-spin bands. The tick marks on the 
energy axes denoted Fe(A,C), Fe(B), V, and Si(p) give the various on-site metal d and metalloid 
p energy levels. Tick marks denoted V in a) refer to the energy levels of a single V impurity in 
FeaX; similarly Fe(B) in b) refers to single Fe(B) impurity levels in Fe2VX. 

FIG. 3. Schematic representation of the common-band model of bonding in a binary AB 
system, a) Densities of states of solids A and B before alloying; b) component densities of states 
in the alloy AB; c) total density of states in the alloy AB. W denotes the width of the common 
band, Wab the bandwidth after bonding. The quantities E^ and E^ give the free-atom energy 
levels while Ea and Eb denote the local on-site energy levels. After Ref. ^. 

FIG. 4. Computed total magnetic moment (per Wigner-Seitz cell) and the moments per atom 
on various inequivalent sites in Fe3_a;Va;Si as a function of V concentration x. All values given in 
units of Bohr magnetons (/U^). Different symbols are explained in the legend. Lines are drawn 
through the theoretical points to guide the eye. Experimental points (open circles) for the total 



moment are after Ref. 14. 



FIG. 5. Same as Fig.0, except that this figure refers to Fe3_2.Va;Ga. Experimental points after 



Ref. 54. 



FIG. 6. Same as FigJ3, except that this figure refers to Fe3_a;Va;Al. Experimental points after 



Refs. 7C and B9 



FIG. 7. Component density of states (CDOS) for various inequivalent sites in FeaSi (topmost 
row). Different angular momentum contributions to the CDOS are shown. The majority (up) and 
minority (down) spin part of the CDOS is given in each case. The dotted vertical lines mark the 
Fermi energy (Ep). Note different scales. 

FIG. 8. Computed partial density of states in the s and d channels associated with a Si site in 
FeaSi; the spectrum has been convoluted with a Gaussian of 2.0 eV (full-width-at-half-maximum) 



to mimic experimental broadening of the soft x-ray emission spectra of Ref. 43 



FIG. 9. Total and site-decomposed density of states in FeaSi, FeaGa, and FeaAl. See caption 
to Fig.^ for other details. 
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FIG. 10. Same as Fig.^ except that this figure refers to Fe2VX alloys. 

FIG. 11. Dispersion curves along high symmetry lines in the Brillouin zone in ferromagnetic 
Fe2VSi (upper) and paramagnetic Fe2VGa (lower). The darker and lighter curves in the upper 
picture represent spin-up and spin-down bands in Fe2VSi respectively. The Fermi energy is marked 
by horizontal lines. 

FIG. 12. Site-dependent density of states in Fes-i^V^Si at Fe(A,C), Fe(B), and V(B) sites for 
each spin direction as a function of the V concentration x. The Fermi energy is marked by dotted 
vertical lines. 

FIG. 13. Total up- and down-spin density of states in Fe3_a;Va;Si, Fes-j^V^Ga, and Fes-ajV^^Al 
for different V concentrations x. The Fermi energies are marked by dotted vertical lines. 

FIG. 14. Same as Figjl^, except that this figure refers to Fea-xV^^Ga. 

FIG. 15. The differences between the total energies of Fes-ajV^^Si and constituent atoms with 
V in the B sublattice (solid curve) and V in the A or C sublattices (dashed curve). 
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TABLES 



TABLE I. The lattice parameters (a) and the muffin-tin radii (Sk) in FesX. 



5fc(A) 



Fes Si 

FesGa 

FesAl 



5.653^ 
5.812^ 
S.TOl'^ 



1.224 
1.258 
1.254 



^Ref. |J. 
^Ref. |2. 
^Ref. 1^. 



TABLE II. Magnetic moments on various inequivalent sites in the end compounds FesX and 
Fe2VX in units of Bohr magnetons (f^s)- The results of the present computations are given together 
with those of a number of other computations and relevant experiments. The values of moments for 
a single V(B) impurity in FesX and a single Fe(B) impurity in Fe2VX are based on our KKR-CPA 
computations in the limiting compounds, and are marked with stars. 





Fe(A,C) 


Fe(B) 


V(B) 


X(D) 


Total 




The 
present 


ory 
others 


Exp. 


The 
present 


ory 
others 


Exp. 


Theory 
present 


The 
present 


»ry 
others 


The 
present 


ory 
others 


Fes Si 


1.34 


1.2^ 
1.10^ 


1.2'' 

1.07<^ 

1.35^ 


2.52 


1.9^ 

2.52^= 


2.4'' 
2.23^^ 
2.20= 


-0.62* 


-0.08 


-0.08^^ 


5.07 


4.63^^ 


FesGa 


1.87 


2.09' 


1.7s 


2.47 


2.38' 


2.5s 


-1.16* 


-0.08 


-0.10' 


6.01 




FesAl 


1.90 


1.36'' 
1.26' 

1.89J 
1.9"^ 


1.5 b 


2.44 


2.48" 
2.16' 
2.34J 
2.4"^ 


2.18^ 


-1.27* 


-0.10 


-O.I4J 


6.03 




FeaVSi 


0.43 


0.37= 




3.08* 






-0.11 


-0.02 


-0.02= 


0.68 


0.70= 


FeaVGa 


0.03 






3.20* 






-0.06 


0.00 




0.00 




FeaVAl 


0.03 






3.18* 






-0.04 


0.00 




0.00 





^Ref. |S. 
''Ref. ||,|5. 
=Ref. H. 
'iRef. 0. 
=Ref. |6. 
^Ref. ||. 
SRef. |2, 
''Ref. |7. 
'Ref. ||. 
JRef. |8. 
'^Ref. Wi. 
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